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Abstract 

We give a tracker solution for the quintessence scalar field for Ratra-Peebles or 
SUGRA potentials, holding before, during and after the equality epoch {pm = Pr) 
and nicely fitting the numerical behavior. 
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1 Introduction 



Dark Energy (DE) is one of the main puzzles of cosmology. It could be a scalar 
field 0, self-interacting through a potential V^(0) (Wetterich 1988, Ratra & 
Peebles 1988), so that 

PDE = Pk,DE + Pp,DE = 0^/2a^ + V{(j)), PdE = Pk,DE - Pp,DE = WpDE 



provided that Pk,DE/y ^ 1/2? so that — 1/3 ^ w > —1. Here 

rfs^ = a^{T){^dT'^ + dxidx') , (i = 1, .., 3) (2) 



is the metrics and dots indicate differentiation with respect to r (conformal 
time). This kind of DE is dubbed dynamical DE (dDE) or quintessence. 

The most significant potentials are those allowing tracker solutions, so limiting 
the impact of initial conditions. Among them, we consider here the SUGRA 
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(Brax & Martin 1999, 2001; Brax, Martin & Riazuelo 2000) and Ratra & 
Peebles (1988, 2003) potentials 

SUGRA : V{(j)) = (A^+"/0") exp(47r0VmJ) (3) 

RP : V{(P) = A^+"/0" (4) 



{nip = Planck mass); in the recent times, they yield |w| close to unity 

and lastly variable or lu^l larther from unity and slowly variable, respectively. 
In a fiat model, once the DE density parameters floE is assigned, either a or 
the energy A can still be freely chosen. 

SUGRA and RP behave differently when approaches rUp. Here, however, we 
shall not deal with late epochs and our expressions, worked out for RP, suit 
SUGRA as well. 

Tracker solutions are then usually found by considering the equation of motion 
4> + 2{a/a)4> = Q;a2A^+"/(/)^+" (5) 



and seeking solutions of the form 

= Ur/nf (6) 



(xj is a reference time and </)j is the field value at that time). They arc fixed once 
the time dependence of the scale factor a(r) is set. In turn, we can however 
write 

a/ a = u{t)It , (7) 



so that -u = 1 (m = 2) in the radiation (matter) dominated eras and is actually 
r dependent around matter-radiation equality. 

In this note we give a single expression for the tracker solution holding before, 
during and after equality, occurring at t^- This expression improves the tracker 

solution expression, known for t <^ t^^ when r approaches Tg . It also neatly 
improves the tracker solution known for r ^ Tg, for a large range of rcdshifts. 
Furthermore, it fits the behavior of the (j) field across equality, with quite a 
small discrepancy from the numerical solution, keeping mostly well below 1 % . 
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2 Tracker solutions through equaUty 



Let he R — a/a^ and 9 — r/rg (og, Tg: equality scale factor, time); it is easy 
to verify that the expression 

R = c^9^ + 2c9 with 0=^2-1, (8) 



gives the scale factor behavior before, during and after equality. This follows 
the integration of the Priemann equation 

Heae{R+iy^^ (9) 



{3He = SnGpe^m = ^T^Gpe,r] Pe,m = Pe,r are uiattcr and radiation energy 
densities at equality) and noticing that H^a^ — 'ic/Te ■ It is then also easy to 
see that 



so that, clearly, m — >• 1 (m — >• 2) for 6' — >• — >• oo). 
Let us now rewrite eq. (5) in the form 
F" + 2^r = ai?2y-(l+a) _ 

9 

Here Y = (^jo with 0"^+° = (aere)^A^"'"" and ' indicates differentiation in 
respect to 9. Let us then seek solutions of eq. (11) of the form 

Y^Y,{R9)'' i.e. (t>^<TY,{aT/a,T,f . (12) 
Clearly 



y = ^ (1 + u)Y , y" = ^ \h\u + If + lm'9 - h{u + 1)1 Y (13) 



so that eq. (11) yields 

Y^+''{R9f^^^°'^ = a{R9fG{9) (14) 

Here 

G-^{9) = + 1)2 + h{u + 1)(2m - 1) + hu'9 (15) 
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Fig. 1. Comparison between a the tracker solution (pir) obtained in this work (GT) 
and the numerical integral (pni'T) of the (/)-field equation. Dotted lines show RT and 
MT solutions vs. the numerical solution. The GT is consistent with the numerical 
integral, apart a narrow interval, where its discepancy from the numerical integral 
marginally approaches 1 % . RT is a fair approximation until Tg. MT performs much 
worse. 

depends on 0, however not so strongly as is meant by a power law. Let then 
be 

6 = 2/(2 + a), Y^^aG, (16) 

although this means that also Y^, depends on 9 and, therefore, the expressions 
(13) should contain further terms. The 9 dependence of G arises from the 
9 dependence of u\ should we neglect the former one, the last term in the 
expression (15) should be consistently omitted. 

In the next section we shall compare the behavior of 0, as obtained from the 
expressions (12), (15), (16), denominated GT, with a numerical solution of 
eq. (5). 



3 High and low 9 limits 

Let us first consider the solution (12) in the limit ^ <^ 1. Then w ~ 1, G ~ 
l/[(26)^ + (26)] and, according to eq. (8), R ~ 2c9. Accordingly, we obtain 

1 

^^'^9^^ , with /3^ = 4/(2 + q;) (17) 



Y. 



Ac^a 



SO that, as expected, we recover the tracker solution for the radiation domi- 
nated regime (herebelow RT). 



a 



Fig. 2. Maximum discrepancy between GT and numerical solutions. 

In quite the same way, for 9 :$> 1, we have w ~ 2, G ~ l/[(36)^ + 3(36)] and, 
according to eq. (8), it! ~ c^O"^. Accordingly, we obtain 



Y 



4 

c a 



1 

2+a 



e^"" , with /3^ = 6/(2 + q;) 



(18) 



and, again, we recover the tracker solution for a matter dominated regime 
(herebelow MT). 



4 Comparison with numerical integrals 



The tracker solution (17) can be used to set the initial conditions to the 
eq. (11), performing then its numerical integration. In doing so we can use 
the expression (8) and (10) so that the system to integrate amounts to 2 first 
order equations. We shall indicate the numerical integral by Yn{6) . In Figure 
1 we plot the ratios Y/Yn, Yr/Yn, Ym/Yn . 

The symmetry between RT and MT solution is only apparent. Dotted curves 
are symmetric in respect to r ~ 10 Tg, and this means that at such redshift, 
inside the matter dominated era, the RT solution still performs as well as the 
MT solution. Using the expressions (17) and (18) it is easy to see that the two 
solutions intersect for 



32 _ 18/3^/2 + 



c2 + 



(19) 



i.e., for = t/t^ ~ vT8/(\/2 — 1) ~ 10 and this confirms what is shown by 
the plot. 



Altogether, this means that a fair tracker solution, since ~ Tg, can only be 
given by using the GT expression, holding until DE itself begins to play a 
significant role as a source of the cosmic expansion. The maximum discrepancy 
between GT and numerical solution is ~ 1 % , slightly after Tg, for any 
reasonable a value. Its a dependence is plotted in Figure 2. 



5 Conclusions 



Tracker solutions play an important role in the analysis of dDE cosmologies. A 
good tracker solution for the radiation dominated regime can be easily given. 
A tracker solution for the matter dominated regime also exists, but its usual 
analytical expression is rather far from a fair numerical behavior. 

Here we gave a fair solution performing as well as the usual one in the radiation 
dominated era, fitting the numerical solution also about the equality (the 

residual discrepancy approaches ~ 1 % just in a narrow interval) and fitting 
the numerical solution all through the matter dominated era. It reads 

= 0e(--)' (20) 



with 



6 = 



2 + q; ' 



(t>e = 



aK{aeTeKf 



h^iu + iy -h{u + i){2u-i) 



(21) 



It was found by using an exact analytical expression for a{t) through the 
equality period, also to work out the r dependence oi u — ra/a . 

When initial conditions to any numerical problem are to be set, one usually 
needs to go back to the radiation dominated era, so to rely on a fair tracking. 
As an example of problems where initial conditions in the matter dominated 
era are useful, let us remind the study of the evolution of a spherical fluctua- 
tion, the prediction of cluster mass function and its redshift evolution (Mainini, 
Maccio & Bonometto 2003, Mainini et al. 2003, Mainini 2005) using a Press & 
Schechter (1974) or similar (Sheth & Tormen 1999, 2002, Jenkins et al, 2001) 
approach, and performing N-body simulations (see, e.g., Klypin et al. 2003, 
Maccio et al. 2004, Solevi et al. 2006). 

Using the GT expression, given here, this can be safely avoided, and initial 
conditions can be given at any time, until DE itself begins to be important 
for the overall cosmic expansion. 
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